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ABSTRACT 

The "Supercritical Pile" is a very economical GRB model that provides for 
the efficient conversion of the energy stored in the protons of a Relativistic Blast 
Wave (RBW) into radiation and at the same time produces - in the prompt 
GRB phase, even in the absence of any particle acceleration - a spectral peak 
at energy ~ 1 MeV. We extend this model to include the evolution of the RBW 
Lorentz factor V and thus follow its spectral and temporal features into the early 
GRB afterglow stage. One of the novel features of the present treatment is the 
inclusion of the feedback of the GRB produced radiation on the evolution of V 
with radius. This feedback and the presence of kinematic and dynamic thresholds 
in the model are sources of potentially very rich time evolution which we have 
began to explore. In particular, one can this way obtain afterglow light curves 
with steep decays followed by the more conventional flatter afterglow slopes, 
while at the same time preserving the desirable features of the model, i.e. the 
well defined relativistic electron source and radiative processes that produce the 
proper peak in the vF v spectra. In this note we present the results of a specific 
set of parameters of this model with emphasis on the multiwavelength prompt 
emission and transition to the early afterglow. 
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Introduction 



The cosmological origin of GRB has by now been firmly establ ished following the 



discovery of their afterg lows and the determination of their redshifts ( I Cost a et al. 



1997 



van Paradijs. et al.lll997l ) and the launch of Swift which increased the number of observed 
afterglows and redshift determinations. These developments left little doubt that GRB 
emission is intima t ely a ssociated with Relativistic Blast Waves (RBW), as proposed by 
Rees fc Meszaros I (119921 ) and at the s ame time shifted the focus of the stu dy from the 
prompt GRB emission to its afterglow (jZhang fc Meszaros 112004 ; |PiranJ|200J). 



The early, sparsely sampled GRB afterglow light curves, were fit well with simple power 



law functi ons, appropriate to emission from either spherical (jSari. Piran. &: Narayan 



19981 ) 



or jet-like (jSari. Piran &: Halpernlll999l ) RBW. However, the launch of Swift with its prompt, 
continuous, broad frequency coverage has provided new unexpected (and unexplained) details 
of the afterglow light curves. Chief amongst them are: (a) An early afterglow steep decrease 
of the flux (oc t~ 3 to t~ 6 ) followed often by a period of constant flux (before its eventual power 
law decline) in many bursts, (b) Large flares in t he X-ray light curves ~ 10 3 — 10 5 sec after the 
beginning of the event (see lO'Brien et al. 1120061 . for more details). These were compounded 
to the already open problems of the prompt emission, namely: (c) The GRB "inner engine" . 
(d) The non-dissipative transport of the GRB energy to the emission region and, most 
importantly, its efficient dissipation there, (e) The physics behind the characteristic energy 
of peak GRB emission, E v and its narrow distribution within the class of the classic GRB 



f lMallozzi et al. lll995tlPreece et al. 1120001 ) . (f) The physics that relate GRB to XRR (X-Ray 



Rich bursts) and XRF (X-Ray Flashes), transients of lower flu x and lower E v , recor ded by 



broad band missions such as BeppoSAX, HETE and Swift (e.g lYonetoku et al.ll2004l ) 



Of the above problems, (a) has received no apparent resolution while (b) is loosely 
attributed to continued activity at the "inner engine"; while not implausible, this demands 
activity over time scales almost 10 7 times longer than the characteristic time associated 
with the "inner engine" dynamics (~ 10~ 3 sec), as the latter is though t related to stellar 
collap se, (d) is considere d to be effected either t hrough protons (e.g iRees fc Meszaros 



19921 ) or magnetic fields (IVlahakis & Konigl I l200ll ). h owever, the necessary and efficient 
dissipation "is one of the least studied aspects of GRB" (IPiran 112004 ): this issue is generally 
approached by parameterizing the energy density in relativistic electrons to be a given 
fraction (typically ~ 50%) of that of protons. Issue (e) is generally open, given the absence 
of an underlying re ason for such a characteris tic energy. Monte Carlo simulations of a large 
number of models (jZhang fc Meszaros 1 120031 ) failed to reproduce the narrow width of the 
observed distribution because of the large number of parameters involved and/or because of 
the lack of strong dependence of E p on any single parameter. Finally, there are a number of 
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prop osals concerning (f) (jYamazaki. Ioka fc Nakamura 1120021 ; iDermer. Chiang, fc Bottcher 
19991 ). which appear plausible but without a single one of them universally agreed upon. 



The "Supercritical Pile" M odel (SPM) (IKazanas. Georganopoulos fc Mastichiadis 



2002; 



Mastichiadis fc Kazanas 1120061 . henceforth KGM02 and MK06), adapted from AGN (IKazanas fc Mastichiad 
19991 ). has been introduced to provide a resolution to (d). The compelling arguments in favor 
of the SPM are: (1) Its economic (non-thermal particles not necessary), efficient conversion 
of the RBW relativistic proton energy into photons through a radiative instability akin to 
that of a supercritical nuclear pile. (2) Its spectra which exhibit a characteristic value for 
E p ~ 1 MeV (in the lab frame) irr espective of the RBW Lorentz factor T, in agreement with 



observation (iMallozzi et al. 1119951 ). produced as "unintended consequence" of the dissipation 
process. Crucial in addressing these issues has been the presence of an upstream medium 
which scatters the RBW photons (a "mirror") and allows them to be re-intercepted by the 
RBW, while boosted in energy by ~ T 2 . 

In MK06 we have explored numerically the SPM assuming a constant Lorentz factor 
T for the RBW, confirmed the efficiency of proton energy conversion into to radiation and 
the presence of a well defined value for E p , reflecting the kinematic threshold of the reaction 
PI ~* pe + e~. The present treatment is far more realistic: (a) It computes the evolution of 
the RBW Lorentz factor T through a medium of density n(r) oc R~ 2 , thought to represent the 
wind of a WR star, including also the effects of the radiative drag of the bulk-Comptonized 
photons, (b) Replaces the upstream "mirror" required by the model by scattering the RBW 
photons in this medium. The combination of these effects can result in a rich GRB time 
evolution, but we presently restrict ourselves to a specific example of a GRB light curve, 
deferring the broader exploration of other models to a future publication. Despite this 
limited scope, we can reproduce some of the salient features of the GRB in the afterglow 
evolution, such as their steep decrease in flux following the termination of their prompt 
phase, an effect traceable in this specific case to the kinematic threshold of the model. 

In §2 we provide the general framework of our model with emphasis on its novel aspects 
compared to previous treatments. In §3 we present the results of our calculations and finally 
in 54 the results are summarized and conclusions are drawn. 



2. The Coupled Radiative — Dynamical Evolution 

We consider a Relativistic Blast Wave (RBW) of speed vq = Pre and Lorentz factor 
T. Its radius R(t) is measured from the center of the original explosion and it is sweeping 
up the CircumStellar Medium (CSM) of density Pcsm- The evolution of T as a function of 
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radius is given by the combination of the conservation laws of mass 



and energy-momentum 



dM n2T1 1 * 

—— = A-kR TpcsM - 
an c 6 T 

dT AirR 2 p CSM T 2 F rad 



(2) 

dR M Mc 2 K 1 

(jChiang fc Dermer Ill999l ). Here E is the radiation emission rate as measured in the comoving 
frame and F ra( j is the radiation drag force exerted on the RBW by any radiation field exterior 
to the flow. Given that the RBW velocity vq is very close to the speed of light c, the entire 
radiative history of the RBW lies just ahead of it at a distance D ~ R/T 2 ; therefore, 
isotropization of this radiation by scattering in the ambient medium (the action of the 
"mirror") will lead to its re-interception by the RBW to thus contribute to -F ra d- This is 
given by the expression 



64vr 



rad 



n nn™ M o*RT*E (3) 
9c 

where r& is the RBW Thomson depth, = Pcsm/^h the CSM electron density and 

<jt the Thomson cross section. In the above expression, two powers of T are due to the 
increase of the photon energy density upon its scattering on the "mirror" while the other 
two to the usual radiative loss rate (an analogous term due to the p7 — > pe + e~ reaction was 
found to increase F rad by 20% for the specific parameter values discussed herein but it maybe 
more important for different values). The calculation of the E and F ra d terms is done by 
implementing the numerical code used in MK06 to compute the radiation of the SPM. This 
is done be solving the simultaneous equations 

^ + L, + Q l = 0. (4) 
at 

The unknown functions n« are the differential number densities of protons, electrons and 
photons while the index i can be any one of the subscripts 'p', 'e' or '7' referring to each 
species. The operators Li denote losses or escape of each species from the system while Qi 
denote injection and source terms of each species by each of a number of processes which are 
described in detail in MK06. The above equations are solved in the fluid frame in a spherical 
volume of radius Rb = R/T. This can be justified by the fact that due to relativistic beaming 
an observer receives the radiation coming mainly from a small section of the RBW of lateral 
width R/T and longitudinal width R/T 2 in the lab but R/T on the comoving frame. 

The present treatment differs from that of MK06 in two important aspects: 

(1) Hot protons accumulate continuously on the RBW as it sweeps the CSM. This 
then sets the source terms of the protons (Q 1 ^) and electrons (Q 1 ^) (with units parti- 
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and 



cles/energy/volume/time) to 

^ nj = £p|(r 2 -r)% P -r) (5) 
QT 3 = PCSM " (r 2 - r)5(ie - r) (6) 

m p m e c 2 R 

i.e. we assume that at each radius R the RBW picks up an equal amount of electrons and 
protons from the circumstellar medium which have, upon injection, energies E p = Tm p c 2 
and E P = Tm P c 2 respective ly. Consequently, the proton energy injection rate is given by 



dBlandford & McKee Ml976h 



^) =47ri?VcsM(r 2 -r)c 3 (7) 

while a fraction m e /m p of the above goes to electrons. 

(2) The scattering of the RBW photons takes place on the CSM ahead of the advancing 
RBW (rather than an ad hoc mirror) and, as such, its photon scattering column is uniquely 
determined by the initial conditions and, like all other parameters is a function of time (or 
equivalently position). 

Eqns. (CE]), ([2]), (@J, along with Eqns. ([3]), ([5]) and form a set which can be solved 
to yield simultaneously the evolution of the RBW dynamics and luminosity. This approach 
is self-consistent in that the 'hot' mass injected through equations ([5]) and shows up at 
RHS of Eqn. ([T]), while the radiated luminosity E feedbacks onto the energy- momentum 
equation through the definition of the radiative force, -F rad , of Eqn. Q. The free parameters 
of this system are (i) the total energy of the explosion E tot (ii) the CSM density profile n(r) 
(iii) the magnetic field as a function of radius B(r). To avoid computation of the evolution 
during the RBW acceleration phase when it likely produces little radiation, we have chosen 
to begin our calculations (and the accumulation of matter by the RBW) at a radius R at 
which is has already achieved its asymptotic Lorentz factor T = T(Rq). 

As proposed in KGM02 and shown explicitly in MK06, the relativistic protons accumu- 
lated in the RBW can become supercritical to the network of jry — > pe + e~, eB — » 7 once 
kinematic and dynamic thresholds are simultaneously fulfilled. The kinematic threshold 
simply reflects the kinematic threshold of the jry —> p e + e~ reaction and reads 

6T 5 >1 or r > 214 (n y 1/12 for B = B eq ~ (8nm p c 2 no) 1/2 T (8) 

with the latter expression for B-field in equipartition, where b = B / B cr a and B cr u = m 2 c 3 /eh 
is the critical magnetic field. The dynamic threshold provides the critical column density 
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Fig. 1. — (a) The Lorentz factor for a RBW propagating in a wind environment with parameters given in 
the text. The thick lines show the evolution with radius without radiative drag, while the thin one with the 
drag included, (b) The multiwavelength spectrum of the burst at times 1, 3, 10, 30, 100 and 300 sec (top 
to bottom). The peak of the bulk Comptonized component is originally close to 1 MeV, however, as the 
evolution proceeds it moves to lower energies. 



for the accumulated relativistic protons to become supercritical (in a fashion analogous to 
a nuclear pile) and, if fulfilled, a large fraction of the energy stored in relativistic protons is 
converted into e + e~-pairs within a few light travel times across the width of the shock. 

At the earliest stages of the RBW evolution the accumulated relativistic proton column 
is small and little emission is possible, only that of the swept-up electrons, which is smaller 
than the energy flux through the shock by a factor m e /m p and may very well represent the 
oft quoted GRB precursor emission. The eventual evolution of the RBW depends on whether 
its asymptotic Lorentz factor To and B— field satisfy the kinematic threshold (Eqn. [8]). If 
not, and in the absence of an accelerated population of particles, only the energy flux in 
electrons is converted to radiation and the GRB is a "dud" , as the combination b T 5 is only 
expected to decrease with radius (however an explosive release is still possible if the proton 
distributi on includes an accelerated power law component that extends to E Tm p c 2 ; as 



hinted in iKazanas. Mastichiadis &: Georganopoulos I (120061 ) , these events may be related to 



P 



the XRRs and XRFs). 

Far more interesting is the case where the kinematical criterion is satisfied initially. 
Then whether the flow becomes radiatively unstable depends on the column of hot protons 
accumulated on the RBW. When this exceeds the critical value, the energy contained in 
relativistic protons is explosively released, the value of -F ra( j increases dramatically and the 
Lorentz factor T of the RBW can decrease over a distance D <C R(t), provided that the 
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Fig. 2. — (a) The bolometric burst lightcurve. The dashed line corresponds to the internal RBW luminosity 
while the dotted line is the bulk Comptonized component. The thick line is the sum of the two. (b) The 
corresponding burst luminosity for various energy bands as a function of time. Long dashed line is at energy 
of 1 MeV, short dashed is at 10 keV, dotted at 100 eV and dot-dashed at 1 eV. The thick full line is the 
bolometric lightcurve. Parameters are as given in the text. 



second term in Eq. ([2]) is dominant. This drop in T is important not only for decreasing the 
emitted flux but, more significantly, for potentially pushing b T 5 below its kinematic threshold 
value ~ 1, (as is the case shown in Fig. CEK), a fact that according to the SPM marks the 
end of the prompt GRB emission phase, i.e. the conversion of proton energy into radiation. 
Following this event, radiation is emitted only from cooling the electrons already present 
within the RBW and those being swept-up by it. The observed flux suffers a precipitous 
decrease with further evolution that depends on the value of R relative to the deceleration 
radius .Rdec corresponding to the resulting value of T; if R < Rd ec , the flux remains at a 
roughly constant level until T resumes its decline, at R > Rdec, through accumulation of 
mass on the expanding RBW; if R > Rdec, then T will continue its decline at the much 
slower conventional level of afterglow theory. 

The detailed, long term evolution of the GRB flux depends on E tot , n(r) and B{r) that 
determine the values of r and V at which the RBW becomes supercritical - it is conceivable 
that for certain parameter combinations supercriticality can be reached at more than one 
radius, with t he released energy being p r oport ional to the time between the corresponding 



bursts see e.g. iRamirez-Ruiz &: Merloni I (120011 ). In Figure [T] and [2] we present the evolution 



of a RBW with n = n (R /R) 2 and B = B (R /R). The parameters are _R = 10 14 cm, 
n = 8.10 s cm" 3 , T = 100, B = 4.4 10 4 G and total isotropic energy E tot = 10 54 erg. 



Figure (T^, depicts the evolution of T as a function of radius in this medium with (thin 
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line) and without (thick line) the radiative feedback. The drop in T corresponds to the 
explosive energy release in the protons and the slow down of the RB W due to the radiation 
drag. As deduced from this figure, R ~ -Rdeo since for R > R , T oc i? -1 / 2 , as expected for 
adiabatic propagation in a wind density profile (thick line). After the decrease in T due to 
the radiative feedback and after the non- adiabatic effects have died out, the evolution of T 
follows a similar track of lower normalization. 

Figure [Tb shows the multiwavelength spectra at various instances as perceived by the 
observer. As it was shown in MK06 the spectrum consists of two components, one that 
is due to the primary particle emission by particles on the RBW and one due to the bulk 
Comptonization of the upstream-reflected primary radiation by the cold pairs of the RBW. 
This latter component peaks early on at lMeV, but as the burst evolves moves to lower 
energies since both T and B drop outward. 

Figure [2^ shows the corresponding apparent isotropic bolometric luminosity as a func- 
tion of time. This consists of the internally produced luminosity (dashed) and that due to 
bulk Comptonization of the mirror-scattered radiation by the RBW (dotted) with the thick 
line representing their sum. As it can also be seen from Fig [Tb, most of the luminosity, 
is by far contained in the bulk-Comptonized component (at E ~ 1 MeV) and exhibits the 
steepest decrease due to the decrease in V and the arrest of additional pair injection from the 
protons. At longer time scales, the only injection available is that of the ambient electrons 
and the emission exhibits the cx t~ l behavior of "standard" afterglows. 

Finally Fig. [2b depicts the luminosity at various energy bands as a function of time - 
here we make no distinction between the direct and the bulk Comptonized component, but 
instead we exhibit their sum. As a rule higher frequencies dominate more at the early stages 
of the burst but drop faster due to a combination of faster cooling and the decrease in T. 
This is consistent with observations: the BAT flux (that receives its major contribution from 
the bulk Comptonized component) decreases much faster than the flux in the other bands 
and its level defines, in effect, the prompt GRB phase (see also next section). 

3. Summary, Discussion 

We have presented above a first attempt at an integrated version of the SPM, complete 
with the coupled RBW dynamics, radiation production and accumulation of hot protons on 
the RBW from the swept-up matter. The latter process is fundamental as the increase of 
the hot proton column to supercritical values is necessary for the explosive energy release 
seen in GRB. Another important feature is the coupling of the radiation to the dynamics of 
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the RBW, the cause of the abrupt decrease in F seen in Fig. [Th. Because this can reduce 
T below the SPM kinematic threshold, it can severely reduce the observed flux, especially 
its bulk-Comptonized spectral component that peaks at E p ~ 1 MeV and constitutes the 
main GRB channel. The existence of the kinematic threshold value for T and its intimate 
association to the radiation emission near E p (~ 1 MeV, the defining GRB property), affords 
for the SPM an operational definition of the GRB prompt phase, a feature unique amongst 
GRB models: as such, the prompt GRB phase is the stage in its evolution during which the 
kinematic threshold condition of Eq. is fulfilled, accompanied by severe reduction in the 
GRB flux following this stage, as observed. 

The time evolution of the flux in Fig. [2] bears great resemblance to that of many Swift- 
XRT GRB, that e xibit a very steep dec lining profile followed by a less steep or flat section in 



their light curves (lO'Brien et al. 1120061 ) . related, as discussed above, to the relation between 
Ro and -Rdec- We believe that the straightforward way that the SPM addresses these vexing 
for the standard model questions attests to its relevance to the GRB underlying physics 
and phenomenology. It should be noted at this point that the efficiency of conversion of 
kinetic energy to radiation depends on the value of the ambient density uq. This dependence 
comes through the dynamic threshold of the SPM, as Uq determines also the value of the 
upstream albedo (i.e. of the "mirror", whose assumption is now obviated), to which the 
dynamic threshold is proportional. We plan to explore the effects of this parameter on the 
GRB properties in a future publication. 

The duration of the burst shown in Fig. [2] is of order of a few seconds. As such it 
would be classified as a short burst, despite the fact that the RBW is assumed to propagate 
in a medium with properties akin to the wind of a WR star. We therefore have presented 
an explicit model that produces a short burst from an object of a young stellar population. 
While it was originally proposed and supported by the earlier observations that short bursts 
are associated with old stellar populat ions (implying neutron star collisions as their source of 



energy), it was shown (IBerger. 1120081 ) that < 1/3 of them are in fact associated with stellar 



populations similar to those of the long GRBs. 

The outlook from this first time- dependent treatment of the SPM replete with the CSM 
distribution and radiation emission and feedback is that this model can potentially produce 
a great variety of GRB light curves (in agreement with GRB phenomenology) which it can 
relate to global parameters of the system. We plan to explore thoroughly these parameters 
in future publications. 
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